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Abstract 

We study the universal fluctuations of the Wigner-Smith time delay for sys- 
tems which exhibit chaotic dynamics in their classical limit. We present a new 
derivation of the semiclassical relation of the quantum time delay to proper- 
ties of the set of trapped periodic orbits in the repeller. As an application, we 
calculate the energy correlator in the crossover regime between preserved and 
fully broken time reversal symmetry. We discuss the range of validity of our 
results and compare them with the predictions of random matrix theories. 
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I. INTRODUCTION 



Over the past decade many studies were devoted to the understanding of quantum man- 
ifestations of classical chaos. This interest can be explained by the fact that this subject 
has applications in many different areas of physics, like properties of complex systems, fun- 
damental aspects of the correspondence principle, transport in ballistic mesoscopic cavities, 
etc.. Most of the theoretical studies have concentrated on spectral properties of closed sys- 
tems, accumulating a large body of numerical evidence of universality and some analytical 
understanding of this fact. Comparatively few studies have so far been devoted to open 
systems, and the scattering problem still lacks some solid fingerprints which serve to clearly 
distinguish integrable from chaotic classical scattering. For this reason, it is desirable to 
study an observable which bridges the well understood quantum aspects of closed chaotic 
systems and the still unclear features of open ones. The Wigner-Smith time delay [|T],0] is 
such an object, since it is intimately related to the level (or resonance) density of the system 
and it is a genuine scattering observable. The present study deals with the universal features 
of the time delay common to all chaotic scattering systems. 

The concept of time delay in quantum scattering was first considered by Eisenbud 
and Wigner |J in the context of one channel scattering. Later on, Smith || extended the 
previous discussions to the many channel problem by introducing the lifetime matrix 

Q ab (E) = -ihY,S ac (E)-^Sl b (E) , (1) 

where S is the standard scattering matrix and the sum runs over all A open channels denoted 
by c. By averaging over the eigenvalues of Q, one arrives at the so-called Wigner-Smith 
time delay 



t(E) = --TV 



S\E)±S ( E) 



-^ l0gdet5(£) , (2) 



which is then interpreted as the typical time spent by the particle in the interaction region. 
Even though this interpretation has some limitations in the case of wave packet scattering 
0], no difficulties arise when the incoming wave can be considered monoenergetic, a common 
situation, e.g. in applications to mesoscopic transport phenomena |||| and microwave cavity 
experiments . 

In general, one can distinguish two regimes associated with a scattering process: a fast 
response (corresponding to direct processes) and a delayed response related to the formation 
of a long-lived resonance. In the energy domain, direct processes rule the energy-averaged 
behavior of t(E). Alternatively, strong fluctuations on the scale of the mean resonance 
spacing A are associated to quasi-bound states, and are, in turn, intimately linked to the 
classical dynamics in the interaction region. 

Our analysis deals with a specific model which illustrates very nicely the most important 
properties of chaotic scattering and is well suited to study the Wigner-Smith time delay. 
Some steps in our considerations take into account system specific properties. However, our 
main results can be easily extended to other chaotic scattering potentials. Our model consists 
of an irregularly shaped cavity (denoted by R in Fig. 1) attached to a pipe (corresponding to 
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the region L) . The boundary between the pipe (or "waveguide" ) and the cavity is arbitrarily 
chosen at the entrance of the cavity, at x = D (the region R need not necessarily be a 
billiard). The quantum propagation in the direction parallel to the pipe axis is free. In the 
transversal direction there are quantized modes 4> c (y) of energy e c , defining the scattering 
channels c. At x — 0, with D chosen to be sufficiently large, the wave function ip{x, y; E) is 
expressed clS 8b superposition of propagating modes 

1>(x, y;E)=J2 (a c e ik ° x - £ S^a^A <j> c {y) , (3) 

c=l V c'=l / 

with the wavenumbers k c given by 

£* = B-u. (4) 

From this equation it becomes evident that by choosing \k\D ^> 1 we ensure that no evanes- 
cent mode survives at x = and Eq. (§) is valid. All the information about the scattering 
process is contained in the energy dependent scattering matrix S(E). 

For chaotic systems, we derive a formula which shows how to calculate the quantum time 
delay using information about the classical orbits trapped inside the system. We further 
explore this formula to compute time delay correlators, which show universal features. The 
universal correlators typically scale with quantities like the average dwell time inside the 
system, Td we ii, and the mean resonance spacing A (also given in terms of the Heisenberg 
time A = 2hTi/ty{). We show that the universal curves obtained in the semiclassical theory 
are in good agreement with the statistical theory of random matrices when r H /rd we ii 3> 1. 

This paper is structured in the following manner. In Section II, a novel derivation of 
the quantum time delay in terms of the underlying classical phase space of the repeller is 
presented. Some applications of this formula are explored in Section III. The comparison 
with random matrix theory is presented in Section IV, where we also discuss the range of 
validity of the correspondence between semiclassical and statistical theories. In Section V 
we present the conclusions of this study. 



II. WIGNER SMITH TIME DELAY AND TRAPPED PERIODIC ORBITS 

This section is devoted to the derivation of a semiclassical equation for the Wigner- 
Smith time delay in terms of classical periodic orbits trapped inside the scattering region. 
Our derivation relies on the association of the S"-matrix with the quantum Poincare map, 
following closely the formalism developed by Bogomolny |10| for closed systems. A similar 



result was previously obtained by Balian and Bloch [JXTj] , based on a construction proposed 
by Friedel |12j for a separable system. The derivation presented below is more transparent 



than the one in Ref. [|Tl|], making the approximations more controllable when dealing with 
actual systems. 

We begin by showing a simple construction relating the energy derivatives of two sets 
of invariants of the A-channel scattering matrix S, namely the eigenphases {8i, 9 2 , ■ ■ ■ , 0\} 
and the traces {Tr S n , n = 1,2, . . .} (the construction remains valid for an arbitrary unitary 
matrix depending on one parameter). For this purpose, let us consider the periodic function 
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F{6) = 9 mod2vr , (5) 



which has a Fourier expansion given by 



F (( ,)=.-2£^. (6) 

n=l n 

Summing over all S'-matrix eigenphases in both sides in Eq.(^) and using ImTrS"™ = 
Y,c=i s ^ n n @c we obtain 

A oo i 

F(6 C ) = Att - 2 Im ]T -Tr S n . (7) 

c=l n=l U 

The convenience of this arbitrary choice of F becomes apparent after differentiating all terms 
in (0) with respect to the energy E, 

*d6 c Ad9 c ^ /n jo n ox ^ldTrS n 
UdE =2K ^dE ^ m ° d ^ ~ 2 Im S ' (8) 

since by recalling (0) it is easy to identify the l.h.s. of the above expression with the Wigner- 
Smith time delay. Eq. (^) was first obtained by Bogomolny |T(J and later rederived by 



Rouvinez and Smilansky ||13|| . Those authors were interested in using the transfer (or scat- 
tering matrix) approach to develop a quantization procedure for closed systems. The ideas 
presented below are quite different, since we are interested in open systems. Indeed, we use 
the closed system to understand the scattering problem, which is the reverse of the procedure 
in Refs. flIB|Jl4]| . To this end, Eq. (R|) only becomes useful after the following steps. 



First, we shall consider the scattering matrix S (at the energy E) for the specific system 
discussed in Section I. Far away from the cavity, at x = 0, the influence of the evanescent 
modes to the wave functions is negligible, since \k c \D 3> 1 (see Fig. 1). Therefore, according 
to Eq. (R|), the exact quantization condition for the system closed at x — 0, becomes det(S' — 



1) = 0, as has already been observed |T4[]. In other words, one of the eigenphases of the 
5*-matrix must vanish 



# c mod27r = 0. (9) 

With this quantization condition, the first term in the r.h.s. of Eq. (Q) is now easily identified 
with the density of states, pl + r(E), of the system closed at x = 0. After proper averaging 
over some energy interval, pl+r(E) can be decomposed into a smooth and a fluctuating part 

A Q/3 

™d2n) = pT +R (E) + pl +R (E) , (10) 

where L and R stand for the pipe and cavity regions respectively (see Fig. 1). 

In the present context, the matrix S is interpreted as the quantum Poincare map of the 
closed system L + R associated with the section S (x = + in Fig. 1). The construction is 
quite obvious, but for sake of completeness let us be explicit: Take an asymptotic incoming 
wavefunction at x = + and let it be scattered by R. As a result, one has a matrix that 
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matches the incoming asymptotic waves into the outgoing ones. This defines the quantum 
return map for x = + and it is also the definition of the S'-matrix. In order to obtain the 
Poincare map, one still needs to describe what happens for x < 0~. This, however, is trivial, 
since x < 0~ defines the asymptotic region (there is no coupling between channels) and the 
corresponding Poincare map is the identity matrix. For the geometry considered here, the 
reflection by the hard wall is equivalent to the reinjection procedure, defining the so called 
Poincare scattering map, originally proposed by Jung [IS|. In this way, the S- matrix can 



also be viewed as the quantization of the Poincare scattering map. It is noteworthy that 
the quantization condition defined by Eq. (|^) has an interesting semiclassical counterpart. 
For closed systems, the accuracy of the semiclassical quantization procedure based on a 
Poincare section requires the section to be traversed by all periodic orbits, such a condition 
defining a "good" section. If this is not the case, evanescent corrections become essential. By 
closing the system sufficiently far away from the cavity region we ensure that the evanescent 
contributions die out, making the exact quantum problem simpler. 

Since the S matrix can be obtained by the quantization of a classically chaotic map, in 
the semiclassical approximation its traces can be expressed as a sum over periodic orbits. 
Actually, the sum over the traces of S in Eq. ([8]) results in an expression very similar to 
the standard Gutzwiller trace formula [TB|] for the oscillatory part of the density of states of 



the system defined by L + R jTJ]. The important difference is that in our case the sum is 
restricted to those periodic trajectories that touch the section S. Decomposing the full set 
of periodic orbits of the system L + R into a set that reaches S and a set which never leaves 
the cavity R, one can write 

1 00 1 f) 

- Im E -^Tr S" * pl +R {E) - p%(E) , (11) 

where p R {E) can be expressed in terms of periodic orbits constricted to the region R. 
Substituting the relations flIU|) and (pTT| ) into Eq. @, we obtain 

^f-r(E) « pf +R (E) + pl +R (E) - (pl +R (E) - p R (E)) , (12) 

yielding 

t(E)k^(pT +r (E)+ P %(E)) . (13) 

This is already very close to the expression we are looking for. The problem is that r is 
measured with respect to x = and we are interested in the time that the particle spends in 
the cavity region, i.e., the time delay with respect to x = D. As a consequence, we still have 
to translate the origin of coordinates to the entrance of the cavity. Under this operation, 
the S'-matrix transforms as 

S(x') = e- ikD S(x)e~ ikD , (14) 

where x' = x — D and k is a diagonal matrix having the fc c 's as elements. Taking into 
account that the time delay is additive with respect to the product of unitary operators, 
namely, 
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Tr((^5 2 )t ^(5x5 2 )) = Tr^-^U) + Tr(4^U) (15) 



we arrive at 



27LD d A, 

T=T "— <16) 

The second term in the r.h.s. of (]T6|) is a smooth function of energy. It is proportional to 
the density of states of the region L (the pipe in Fig. 1), 

2HD d * 2D,£ 1 2nh 

— dBS*- = xS^ = ir^- (17) 

The first identity above says that under a spatial translation by D, the time delay varies 
by the classical time it takes to travel to the point displaced by D and back, averaged over 
the channels. The second identity is obtained by using the Weyl formula for the quasi one- 
dimensional density of states in a waveguide, expressed in terms of the longitudinal velocities 
v c = hk c /m. Inserting ([R]) and flI7D into (0), we obtain 

2lTh/ p%(E)+ P l(E)) . (18) 



A 

This is the main result of this section. The first term in the r.h.s. of Eq. (|18"D represents 
the mean time spent in the cavity, (r) = 27r7i/(AA) = t#/A, in agreement with Levinson's 
theorem (which holds irrespective of whether the underlying classical dynamics is chaotic or 
not; see e.g. [fL7|1 ). The second term is given by |j 



r fl {E) « - Y,Tu{E)A um {E)e imSv{E)/h -^^ m , (19) 

urn 

where the sum runs over the primitive periodic orbits v which do not leave the cavity and 
their m-th repetitions. As usual, one needs from each periodic orbit v its period T v , action 
s, Maslov index fi u , and the amplitudes A v given in terms of the monodromy matrix M v 

B 

A V m = — ; • (20) 

v /|det(M--l)| 

We base all considerations that follow on Eq. (|P9"D . One has to keep in mind that Eq.fliDp is 
a semiclassical result, sharing most of the usual limitations of the standard trace formula for 
bound systems. It must be also emphasized that for the quantization of the bound systems, 
one is free to seek a Poincare map which includes all the periodic orbits (in our system this is 
the Birkhoff bounce map). Then, the r.h.s. of Eq. (|TTD is zero and the "delay time" for this 
map is just the smooth Weyl density of states as discussed by Rouvinez and Smilansky ||T3f . 
The scattering map is predetermined in our case. This is a chaotic, singular, discontinuous 
classical map because of its exclusion of the repeller, which accounts for the fluctuations of 
the time delay. For this reason we cannot use the time delay calculated directly from the 
semiclassical scattering map, reworking the theory on the basis of the repeller. 
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III. UNIVERSAL CORRELATIONS IN THE QUANTUM TIME DELAY 



Here we apply the results of the last section to analyze the fluctuations of the time delay 
as a function of energy. We study the crossover regime in which time reversal symmetry is 
lost due to the presence of an external magnetic field B. Our derivation extends previous 
results by Bohigas and collaborators |T8j (closed systems at the crossover regime) and by 
Eckhardt |19| (open systems with preserved time- reversal symmetry). Our presentation 
closely follows the discussion in [TBI . 



A usual measure to characterize the time delay fluctuations is the correlation function 
K T {e,B) = (rf l (E + s,By i (E,B))^ , (21) 



where (...)# stands for energy average. The semiclassical approach to the calculation of this 
correlator begins by inserting ([19]) in (|2T|) (the dependence of all quantities on B will now 
remain implicit to simplify the notation): 

K M = ^ ^ (T,(E + e)T I/f (E)A I/m (E + e)At lm ,(E) (22) 



v m u'm' 



x exp [ims u (E + e)/% — m[i vm j1 — im! s v i{E) /h — m{i v i m i /2] J E . (23) 

Here the average is taken over an energy interval 5E which, to be meaningful, must be 
large as compared with the quantum scale A in order to include many resonances. On the 
other hand, for practical purposes, SE must be small enough to allow for the use of classical 
perturbation theory. According to these considerations, the most important effect of varying 
the energy is on the actions, as they are measured in units of h. We write 

s(E + e) « s(E) + eT , A(E + e) « A(E) , and T(E + e) « T(E) . (24) 

Next, to evaluate K T (e), we use the "diagonal approximation", neglecting contributions 
of pairs of orbits having distinct actions, as they cancel out upon averaging over energy. This 
approximation is accurate for orbits with periods shorter than some critical value, whose 
typical scale is the Heisenberg time th- We shall come back to this point later. 

For B = only two kinds of pairs of orbits survive: the pairs of identical orbits and the 
pairs of time-reversed partners. As the magnetic field B grows, the contribution to the cor- 
relation function of time-reversed partners gradually decreases. Keeping both contributions, 
we have 

= 4 E ( T u \A um \ 2 e iemT "' h [l + e im5s ^ h ] ) e , (25) 

um 

where 5s = 5s(E, B) stands for the action difference between a pair of time-reversed orbits. 
Next, we group orbits having periods in a small interval 5t around t (containing many 
orbits). This introduces a kind of averaging procedure, defining smooth functions of t. We 
then integrate over all t 

K ^=^JZ dtltlei£t/h ( t|i42 )*([ i+e<WR ])t ) ' (26) 
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where we discarded the multiple repetitions \m\ 7^ 1, which are exponentially negligible with 
respect to primitive orbits. 

To evaluate the time average of the amplitudes A we use the sum rule for open chaotic 
systems [p0|j21]l : 

)t\A 2 ) = e -?™ , (27) 



where j(E), the so called escape rate, is 7 = l/rd we ii, with rd wc ii defined in the Section 
I. Eq. ( p7|) has a simple physical interpretation: As we increase the period, periodic orbits 
proliferate exponentially (with a rate given by the entropy) as their stability tend to decrease 
also exponentially (the rate given by the sum of positive Lyapunov exponents). For 7 = 0, 
both effects cancel and one recovers the sum rule based on the uniformity principle for closed 



systems [p2 |. For open systems the disbalance between the entropy and Lyapunov exponents 



is just the escape rate resulting in ( |27| ) p3| , [2"l 



The time average of the crossover term has been discussed in detail in Ref. ||18|| . Further 



theoretical arguments []25| and numerical evidence [^6| suggest the exponential decay 

( 1 + e ^/n^ = 1+e -B^(E)\t\/K* (2g) 

Here k(E) is a purely classical quantity, which measures the rate of decay of the appropriate 
correlations in the chaotic system. For more details see Ref. jnj. Identifying 'y(E) and k(E) 
with their mean values (the energy averaging interval is small) we write 

K T {e) = r dt \t\ e l£t ' h -^ ( 1 + e -B 2 <E)\t\/^\ (2Q) 



This integral is easily evaluated, resulting in 

*» = 2< r > {FTW + [(r + ^ + ^ J ' (30) 

where, for the sake of future comparisons, we have defined 

lu = ite/A, r = 7r 7 n/A, y=?J?L. (31) 



When B = 0, we note that Eq. (|30|) reduces to the result obtained by Eckhardt |19| for the 
time reversal symmetric case. Alternatively, putting 7 = we recover the results of Ref. 
18| for the density-density correlator of the closed problem. 

One of the most interesting properties of the semiclassical approximation to K T (cu) is 
that in this case we obtain an intrinsically more accurate result than the density- density 
correlators for closed systems, that have been extensively semiclassically studied. The reason 
is that the semiclassical approximation starts to fail for energy domains of the order of the 
mean level spacing A, corresponding to times longer than th = 2-n-h/A, and such times are 
normally unimportant for the computation of K T . The physics of the scattering problem 
provides us with a natural cut-off for the summation over periodic orbits, which is given by 
the typical escape time I/7, explicit in (|27J). Since the semiclassical theory is only applicable 
if the waveguide L has many open channels, I/7 must be much smaller than the Heisenberg 
time th, corresponding to the regime of strongly overlapping resonances. Although we do 
not have a control over the magnitude of the accuracy, by increasing the number of open 
channels in actual systems, K t (uj) converges to the semiclassical approximation 
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IV. COMPARISON WITH THE STATISTICAL APPROACH 



The statistical properties of the Wigner-Smith time delay for chaotic systems have been 
also investigated by using the Random Matrix Theory |f2~T|-|31"|| . In analogy with the semi- 
classical approach, the statistical approach requires a decomposition of the Hilbert space 
into an asymptotic region L and a "complex" scattering region R, following the notation of 
Fig. 1. As before, depending on the energy E, one will have A propagating channels in the 
pipe, labelled by c. By introducing arbitrary boundary conditions at x = D one can define 
a set of quasi-bound states \fi) in region R and a set of scattering states \x c (E)) in L. These 
states form a complete set Q + P = 1, with Q = J2fi l/^X/^l an d P = J2 C J dE\x c (E)) (x c (E)\. 
Thus, the Hamiltonian is given by 



H = 



QHQ + PHP + QHP + PHQ = H QQ + H PP + H PQ + H QP , (32) 

where Hqq is interpreted as an "internal" Hamiltonian, and H P q (Hq P ) are the couplings 
between the interior R and the channel region L . After some algebra, one can show |32|JT7|| 
that the resonant S-matrix can be written as 

S = I - 2mH PQ 1 H QP , (33) 

& - n QQ + VKtlQp n PQ 

in the absence of direct reactions, implying that H PP is diagonal. The decomposition of the 
Hilbert space in projectors P and Q can, in principle, be employed for a large variety of 
problems, which makes Eq. fl3"3| ) a very useful parameterization of the S'-matrix. It follows 
that the Wigner-Smith time delay is given by |3l| 



r{E) = ~ImTr(E - H QQ + mH QP H PQ ) \ (34) 

This expression is akin to the one obtained semiclassically QI%D, as it should. Here r(E) 
is equated to the level density of the "closed" system, defined by the operator Q (which is 
quite arbitrary) and smoothed by the coupling to the exterior world by the imaginary term 
in (0). Although conceptually similar to Eq. (|i~8|), it is not a simple task to arrive at the 
semiclassical expression starting from Eq. (|34l) . 

Since one expects signatures of chaos in scattering processes to be manifest for times 
much longer than the typical traversal time, we focus our attention only in the resonant 
part of S. This is the physical justification for neglecting direct (fast) reactions in Eq. (P3"|). 
Moreover, the object which is responsible for classical chaos in scattering is the repeller, 
implying that chaos is a property of the "internal" Hamiltonian. Therefore, in analogy with 
Bohigas' conjecture for closed systems, a statistical modelling of quantum chaos in open 



systems can be made by taking Hqq as a member of an ensemble of random matrices |L7 
For instance, for preserved time-reversal symmetry Hqq belongs to the Gaussian Orthogonal 
Ensemble (GOE) and for broken time-reversal symmetry to the Gaussian Unitary Ensemble 
(GUE). This conjecture allows us to study universal fluctuations in scattering processes 
by calculating S'-matrix correlation functions. Those are obtained by ensemble averaging, 
which is equivalent to an energy averaging based on the ergodic hypothesis. 

In particular, the calculation of the 2-point time delay correlation function K T (e,Y), 
studied in the previous section, requires the evaluation of 
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K T (e, Y) = A Re | Tr g(E + |, Y)Tr g^E - |, Y) - Tr g(E, F) 2 } , (35) 

where O denotes the ensemble average of O and 

g(E,Y) = (E — Hqq{Y) + mHqpHpo) 1 , (36) 



with the variable Y parameterizing changes in the internal Hamiltonian. For instance, if 
Y stands for an external magnetic field B, one can study K T (e,Y) in the crossover regime 
between preserved and broken time reversal symmetry by choosing Hqq = H GOE + Y H GUE . 
The results are universal in terms of the scaled variable y = Y/Y c , where Y c is system specific. 

From the technical point of view, the ensemble average in Eq. (|3"5] ) implies a nontrivial 
calculation based on the supersymmetric technique developed by Efetov [[34[] . This technique 



was adapted to scattering problems by the Heidelberg group p5| , p6| . Ref. |3oj is the starting 
point of all works that use the statistical approach to study the time delay correlation 
function |3(],[37|,[n|] and related objects f38|. A discussion of the supersymmetric technique 



is beyond the scope of this paper and we refer the reader to the excellent introductory text 
in Ref. f39j, and to the S'-matrix review in Ref. pTj . 



Let us start with the simplest case, K T (uj,y ^> 1) for broken time-reversal symmetry, 
corresponding to taking Hqq as a member of the Gaussian Unitary Ensemble. The result, 
given as usual in terms of a double integral, can be found in Ref. |31| 



*f» = { 4- f\ dx f ^ c - - A )) n ( 2+^-^ ) ^ 



1 Ji 



where the transmission coefficient T c gives the probability of an incoming wave at the channel 
c in the vicinity of x = D to enter the scattering region R. In order to compare with the 
semiclassical theory, one has to take T c — 1 for all channels c, since this theory does not 
take into account any barriers preventing perfect transmission. Thus, keeping the notation 
introduced in Section III and identifying (r) = 1/7, the leading asymptotic term in powers 
of T _1 of Eq. ( PV| ) becomes 

K gue (uj) w f2 (w) 

in nice agreement with the semiclassical result. In Fig. 2 we present K GUE as a function of 
the number of open channels A. As A increases, the agreement with the semiclassical theory 
becomes much better, as it is nicely shown in the inset of Fig. 2. Even for relatively small 
A, the exact result does not differ significantly from the semiclassical one. This is explained 
by the fact that one can show that the next to leading order correction is smaller by a factor 

r 2 . 

The other simple limit is the case where time reversal symmetry is present, corresponding 



to Hqq taken as a member of the GOE. Here the result is |30| 

(r) 2 * 



(t) r 1 f°° f°° 
K goe (uj) = — ^— dX d\ x rfA 2 /i(A,A 1 ,A 2 )(2A + A 1 + A 2 ) 2 (39) 
4 jo Jo Jo 

(1-T C A) 2 
l + T c Xi)(l+T c X 2 ) 



cos ( W( 2A + a, + a 2) ) n f n ^rJ;T T ,.J 1/2 . (40) 
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with the measure fi given by 



, (i_A)A|Ax-A 2 | 

M ' Al ' 2j [(1 + Ai)Ai(l + A 2 )A 2 ]V2( A + Ai) 2 (A + A 2 ) 2 ' 1 j 

In this case, even taking T c — 1 for all c's the integral is still difficult. However, we can use a 



trick introduced by Efetov [gj] or the asymptotic expansion proposed by Verbaarschot |36 
to obtain 

again in nice agreement with the semiclassical result. Here, higher order corrections are 
smaller only by a factor T. This is manifest in Fig. 3, where one can see that the GUE case 
converges faster than the GOE to the semiclassical result. 

For the crossover regime the supersymmetric expressions become even more complicated. 
By performing a calculation similar to the one done by Pluhaf and collaborators [f40| |, Fy- 



odorov, Savin, and Sommers j37| obtained a closed expression for K T , first numerically 



studied in Ref. |41| as the ballistic limit of electronic mesoscopic transmission. The leading 



asymptotic term in inverse powers of T is 



KMv) = if jp^p + [ ( r4^p ) • («) 



identical to the semiclassical result Eq. 

The strength of the statistical method is in dealing with situations beyond the scope of the 
semiclassical theory, either by analyzing situations where A is small, or by treating systems 
where the scattering waves have to overcome barriers. In the latter case, the semiclassical 
analysis needs some refinements. Equally or even more problematic is the fact that when 
J2T C < 1, the resonances become isolated. This calls for a semiclassical theory that deals 
with times larger than r#. 



This is a good point to call the attention to one issue in the literature |23] that is often 
confusing. The ensemble average is performed for a fixed number of channels and many 
resonances between channel thresholds. The accuracy of the statistical approach increases 
as the number of resonances increase. With such a construction it is not possible to use 
this theory to statistically analyze the S"-matrix of a system like the 3-disc problem, since 
between channel thresholds there are very few resonances, explaining the disagreement in 
|24| . However, since one takes the trace over the channels to obtain the Wigner-Smith time 



delay, provided A>l,a small change in A seems not to affect K T , as observed in Ref. [ITU . 



V. CONCLUSIONS 

In this paper we presented a semiclassical derivation of the formula connecting the 
Wigner-Smith time delay r to the resonance density of the scattering region, corresponding 
to a chaotic system. We showed that r can be written as a sum over the periodic orbits 
inside the repeller. The physical interpretation of this relation is that the repeller is respon- 
sible for the time spent in the cavity by the scattering trajectories. An open trajectory that 
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closely approximates a periodic one, can spend a long time in the scattering region. This 
dwell time essentially depends on the stability of the periodic trajectory that is approached. 
As a result, the typical classical dwell time depends on few bulk characteristics of the scat- 
tering system. The interesting achievement of the semiclassical theory is that it is possible 
to write scattering observables in terms of classical trajectories that never leave the system. 
Doing so, one avoids all problems inherent of a semiclassical formulation in terms of open 
trajectories. 

One of the striking features of the semiclassical approximation for a scattering system 
is its accuracy. In distinction to the usual studies of density- density correlators in closed 
systems, here one has no need to account for times of the order of r H , the time scale where the 
semiclassical approach starts failing. Trajectories entering an open chaotic system cannot 
typically stay inside the scattering region for times much longer than (r). This fact provides 
us with a natural cut-off for any semiclassical summation formula, namely (r) itself. If 
we assume that convergence is the only problem of the semiclassical formalism, systems 
with increasing numbers of open channels will be described by the semiclassical theory with 
increasing precision as compared with the exact theory. 

Although we already know the exact statistical result for several correlators and distri- 
butions involving the Wigner-Smith time delay, such an approach, by construction, does 
not have information about non-universal quantities (like A, Y c , etc.). Those are usually 
extracted from the experiment. The point of this paper is that this information is usually 
available from the classical dynamics and the semiclassical approach can always be adapted 
to give a recipe to compute the non-universal scaling factors. In summary, even if the semi- 
classical theory cannot compete in accuracy, it can be used as a complement to the statistical 
approach. 
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FIGURES 



FIG. 1. Schematic illustration of the scattering system under investigation. The "cavity" region 
is denoted by R and the attached waveguide by L. 

FIG. 2. Comparison between the semiclassical (solid lines) and the exact random matrix results 
for the time-delay correlator (case where time-reversal symmetry is absent). We plot the normalized 
correlators K* = 2T 2 K T / (r) 2 vs. the normalized energy co/T. Different line styles indicate different 
number of open channels; dashed, dash-dot, and dotted correspond to A = 3, 5, 10, respectively. 
Inset: difference between random matrix results and semiclassics. 

FIG. 3. Comparison between the semiclassical (solid lines) and the exact random matrix results 
for the time-delay correlator (case of preserved time-reversal symmetry) . We plot the normalized 
correlators K* = T 2 K t /(t) 2 vs. the normalized energy uj/T. Different line styles indicate different 
number of open channels; dashed, dash-dot, and dotted correspond to A = 3, 5, 10, respectively. 
Inset: difference between random matrix results and semiclassics. 
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